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^ Abstract 

^ We use an effective field theory for short-range forces (SREFT) to analyze systems of three 

Q identical bosons interacting via a two-body potential that generates a scattering length, a, which 

is large compared to the range of the interaction, £. The amplitude for the scattering of one boson 
off a bound state of the other two is computed to next-to-next-to-leading order (N^LO) in the 

O i/a expansion. At this order, two pieces of three-body data are required as input in order to 

renormalize the amplitude (for fixed a). We apply our results to a model system of three Helium-4 
atoms, which are assumed to interact via the TTY potential. We generate N^LO predictions for 

> 

atom-dimer scattering below the dimer breakup threshold using the bound-state energy of the 
00 shallow Helium-4 trimer and the atom-dimer scattering length as our two pieces of three-body 

input. Based on the convergence pattern of the SREFT expansion, as well as differences in the 
^Sj predictions of two renormalization schemes, we conclude that our N LO phase- shift predictions 

will receive higher-order corrections of < 0.2%. In contrast, the prediction of SREFT for the 

• ^ 

^ binding energy of the "deep" trimer of Helium-4 atoms displays poor convergence. 



Dedicated to Professor Henryk Witala at the occasion of his 60th birthday. 
Permanent address 



Electronic address: jichen@triunif.ca 
Electronic address 



phiilipsOphy.ohiou . edu 



I. INTRODUCTION 



Few-body systems share universal features at low energies, for which details of their short- 
distance structure and interactions are not essential. Studies of such systems' universal 
behavior normally involve separation of a large-distance scale and a short-distance scale. 
The large-distance scale is related to the two-body scattering length, a, which determines 
the zero-energy total cross section in S-wave elastic scattering of two particles: 

a = Ana^ , (1) 

where a = 27^0? if the two particles are identical. The short-distance scale is represented 
by the range of the interparticle short-distance interactions, £. Alternatively, the separation 
of scales is also reflected in a hierarchy of momentum scales. The low-momentum scale, 
Q ~ 1/a, is the typical momentum scale of such few-body systems (e.g. the binding mo- 
mentum). A high- momentum scale, ~ sets the breakdown scale for the description of 
universal physics: high- momentum degrees of freedom above l/£ (~ short-distance physics 
corresponding to interactions of range are integrated out. 

Universal behavior thus occurs in few-body systems where the two-body scattering length 
is much larger than the range of the interaction, \a\ ^ ^. For example, in systems with three 
equal-mass particles, Efimov proved the existence of an infinite tower of three-body bound 
states (trimers) in the unitary limit \a\ — )■ oo [1]. These trimer states have a geometric 
spectrum: the ratio between two consecutive binding energies is always 515. This behavior, 
along with many other universal features, has been studied in systems in atomic, nuclear and 
particle physics which obey |a| ^ £. In ultracold atomic gases, the atom-atom scattering 
length is controlled by an external magnetic field and can be tuned to arbitrarily large values 
through Feshbach resonances [2j. A large scattering length occurs without manipulation in 
the interactions of Helium-4 atoms. Experimental data [3] suggests that the ratio a/£ is 
of order 10 in this system. In few-nucleon systems, the nucleon-nucleon (e.g. np and nn) 
scattering length is naturally about 3 times the range, a/£ ~ 3 These systems, together 
with halo nuclei [SI E] and exotic charmonium [7j, offer testing grounds for the study of 
universal physics. For reviews of the application of this notion of universality in few-body 
physics, see Refs. [HI E] . 

One powerful tool to describe universal physics is Effective Field Theory (EFT). It pro- 
vides a systematic expansion in the ratio of a low- and a high-momentum scale. It thereby 
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permits study of not only the universal features at low energies, but also corrections to these 
features from dynamics at (or below) the length scale which can be evaluated in pertur- 
bation theory. For example, the binding of nucleon-nucleon systems is much smaller than 
the pion mass, m^. Such systems' low-energy behavior is thus insensitive to physics with a 
momentum scale at or above 771^^, e.g. pion exchanges. It can be described by an effective 
theory, called pionless EFT, in which nucleons are treated as point-like particles and the 
nucleon-nucleon potential becomes a string of contact interactions with increasing powers 
of momentum. These contact terms parametrize short-distance physics in nucleon-nucleon 
systems, and all those with two or more derivatives are investigated perturbatively when 
higher-order corrections are examined in this EFT. 

Pionless EFT has been successfully applied to nucleon-nucleon systems and is able to 
describe the low-energy component of nuclear forces [T0Hl5]. Using renormalization-group 
methods, Birse et al. showed that, if \a\ ^ i, this EFT is equivalent to the effective- 
range expansion for scattering by a short-ranged potential (TB]. Therefore, this EFT can be 
understood as an expansion in the ratio of the effective range, tq, to the scattering length, 
a. These are defined from the effective-range expansion of the two-body S-wave phase shift: 

kcot6 = + -rok'^ + . . . . (2) 
a 2 

We refer hereafter to this EFT as short-range EFT (SREFT). If spin and isospin degrees 
of freedom are excluded, SREFT can be used to describe low-energy boson-boson systems 
with a short-range interaction. 

The SREFT expansion is valid for k ~ ^ 1/^, and tq is usually of the same order as 
i. When \a\ ^ £, leading-order (LO) calculations in SREFT are computed in the zero-range 
limit i = 0, and a is treated as an input quantity fixed by extrapolation from scattering 
data or by model calculations with an underlying interaction beyond EFT. Effects of the 
short-distance physics are systematically included in SREFT as higher-order corrections to 
the zero-range limit. At these orders tq is also an input, to be obtained in the same way as a. 
Once a and tq are fixed, results from these EFT calculations are insensitive to short-distance 
details. 

The universal features of three-body systems, which were predicted by Efimov, have 
recently been rederived in SREFT by Bedaque et al. [I7l[THj. The leading-order calculation 
corresponds to the limit i = and large a. In SREFT at LO, diagrams involving loops still 
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contribute to the same EFT order in powers of momentum as do tree diagrams, which thereby 
leads to a non-perturbative calculation for three-body observables in this EFT at LO. A 
three-body counterterm is needed at LO for consistent renormalization, i.e. to cancel notable 
dependence of the non-perturbative result on the momentum cutoff, A. This counterterm 
is introduced to describe short-distance physics in the three-body system, and is tuned to 
reproduce one three-body observable. 

Beyond LO in SREFT, three-body observables are calculated as a perturbative expansion 
in powers of r^/a. The next-to-leading-order (NLO) range corrections, ~ rg/a, have been 
calculated in [121 ED] for systems with a fixed scattering length, and in [2TH23] for systems 
with a variable scattering length. The latter is associated with needs in the experimental 
study of cold atoms, where the atom-atom scattering length is manipulated near a Feshbach 
resonance. The existence of a scattering-length-dependent parameter in the NLO three-body 
counterterm was reported in [22l [23] . 

Meanwhile, few-body systems with a large, fixed scattering length are realized in nuclear 
and molecular physics. In the nuclear case r^/a ~ 1/3, and so we expect NLO predictions 
to have errors ~ 10%, which is above the desired accuracy when comparing to experiments. 
(For example, even the thirty-year-old measurement of the neutron-deutron scattering length 
by Dilg et al. has an accuracy of about 6% [21] •) Therefore, next-to-next-to-leading-order 
(N^LO) effective-range effects, ~ need to be considered in nuclear systems. 

Such an N^LO calculation was carried out by Bedaque et al. [20], who calculated the 
S-wave neutron-deuteron phase shift at N^LO and showed that an additional, energy- 
dependent, three-body counterterm is needed in the renormalization. This conclusion is 
supported by analytic arguments based on a renormalization-group treatment of the three- 
body problem in hyperspherical co-ordinates [25j. However, the calculation of Bedaque et 
al. is not done according to a strict expansion in ro/a. Instead, range corrections up to 
0{rQ/a^) are included in the two-body scattering amplitude, which is then used as an in- 
put in a non-perturbative three-body calculation, thereby arbitrarily including higher-order 
corrections above N^LO. Bedaque et al. argue that these higher-order corrections should be 
small, as long as the cutoff A is kept below l/£. 

Platter and Phillips analyzed the three-boson system up to N^LO by using a similar 
partial resummation as in Ref. [20], and concluded that an energy-dependent three-body 
counterterm is not needed for renormalization in the limit A ^ l/rg [26j. This contradicts 
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Bedaque et a/.'s result [20]. However, it must be emphasized that Ref. [26] only reached 
this conclusion for the case that the cutoff A — )• oo. Ji et al. [23j subsequently showed that 
the partial resummation of NLO range corrections carried out in Ref. PS] is only consistent 
with a rigorous perturbative expansion when I ^ tq, with A ~ l/l. This condition should 
also apply to calculations at N^LO. 

In order to solve the controversy raised by these works, a N^LO calculation in a rigorous 
perturbative expansion of three-body observables is needed. In this paper, we calculate the 
Tq/o? corrections to three-body bound-state energies and scattering phase shifts perturba- 
tively, and compare our findings with those two previous works [201 EH]- Our derivation is an 
extension of the NLO calculation reported in [23]. In the rigorous perturbative expansion, 
we find that an additional energy-dependent three-body counterterm needs to be included 
for a consistent renormalization at N^LO. This result thus supports Bedaque et al.'s con- 
clusion, and shows that Platter and Phillips' findings do not apply to the generic situation 
where tq ~ 

In Sec. [TT], we will briefly review the SREFT in a three-boson system, and summarize the 
formalism at LO in [17J, which is based on a modified Skorniakov-Ter-Martirosian integral 
equation, as well as the formalism at NLO, which is based on a perturbative expansion [T9l 



. We then discuss in Sec. Ill the calculation of the N LO three-body scattering t-matrix 



and relate it to the real K-matrix amplitude. Following that, in Sec. IV, we relate that 
K-matrix to an N^LO calculation of three-body observables: three-body binding energies 
and the S-wave phase shift. Sec. |V] shows the asymptotic behavior of the LO three-body 
amplitude and defines several relevant integrals. This section serves as background for 



Sec. VI, which discusses renormalization at N^LO, via an analysis of the regularized results' 
cutoff dependence and the counterterms needed to cancel it. We also show our numerical 
result for the N^LO three-body force as a function of cutoff and energy, and compare it 



with analytic expressions. In Sec. VIL we apply our formalism to observables in Helium 
trimers, comparing our results with precise calculations done with both SREFT and model 
potentials in Refs. [261428] . We conclude with a summary, followed by appendices. 
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II. THE SREFT IN THREE-BODY SYSTEMS 



In this paper, we simplify our consideration to a non-relativistic system including three 
identical bosons. In particular, the simplifications in bosonic cases do not result in a loss of 
generality, due to the universal behavior in non-relativistic three-body systems. As shown 
by Bedaque et ai, the formalism for three identical bosons can be straightforwardly adapted 
(with appropriate modifications for the inclusion of spin and isopsin) to the study of three- 
nucleon systems ^j. Therefore our N^LO study in bosonic systems can be extended to 
three-nucleon systems, such as the neutron-deuteron system (including ^H) and the proton- 
deuteron system (including '^He). (Although, in the latter case, the Coulomb potential 
between two protons must be accounted for [301 SI]-) For terminological simplicity in our 
three-identical-boson system, we from now on refer to a single boson as a (bosonic) atom, a 
two-boson state as a dimer and a three-boson state as a trimer. 

Because we are considering only short-range interactions, we write the Lagrangian in the 
SREFT to describe non-relativistic three-atom systems, so it includes an atom field (ip), a 
dimer field (T), an atom-atom contact interaction and atom-dimer contact interaction: 

£ = ^t (id^^^—\ ^ + aT^ (idQ + — - /\ \ T--^ (TV^ + h.c)+/iTtrV^V + --- • (3) 
\ 2m / \ Am / v 2 

The ellipses represent interactions above LO in powers of momentum, that are suppressed at 
low momenta. As demonstrated by Kaplan [32] , a positive effective range can be described 
by the theory with a = — 1, which describes an atom and an auxiliary dimer. 

In the following, we will employ Eq. ^ and expand each quantity described by this 
Lagrangian in powers of tq. In the tq expansion, this Lagrangian is order-by-order equivalent 
to one containing only a single atom field (for detailed proofs of this equivalence, see Refs. [T71 
[T8l [33]). The atom-dimer contact interactions in our SREFT are equivalent to the three- 
atom counterterm in Refs. [171 HE] , which is an equivalent description of the short- distance 
physics in a three-body system. Thus, from now on, we will refer to the atom-dimer contact 
interaction as a three-body counterterm (or three-body force). 

One might be concerned that the van der Waals' potential which governs the longest- 
range part of the atom-atom interaction is not truly short-ranged. It is indeed the case 
that the effective-range expansion for atom-atom scattering via a van der Waals' potential 
contains a term ~ fc^, i.e. a term that is non-analytic in fc^. Gao [3^] showed that, upon 
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including the van der Waals interaction, the two-body S-wave effective-range expansion is 
modified to: 

kcot6o = -- + JroP - Prlk' + O f"^] . (4) 
a 2 \ J 

However, the term is also suppressed by r^/a relative to the shape-parameter term ~ fc"^, 
and so is less important than that piece of the effective-range expansion in either of the 
regimes ~ 1/a or A; < l/'^o- Since the shape-parameter term does not enter SREFT 
calculations until N^LO we conclude that SREFT can describe atom-atom interactions up 
to at least N^LO. 

The atom propagator is formulated as 

^'5(po,p) = -2 r> (5) 

where po and p denote, respectively, the time component and the space component of 4- 
momentum. The dressed dimer propagator is renormalized to reproduce tq and a in the 
effective-range expansion for the atom-atom scattering amplitude (see Eq. (|2|), and is ex- 
pressed as 

zP(po,p) = ] , , (6) 

-7 + 2^0(7^ + - p^/4) + ^J-nipQ + p2/4 - ie + ie 

where 7 is the "typical momentum" of the atom-atom system, which is, from Eq. (|2|, related 
to a and tq by: 

i.,_i.„,. + o(r|). (7) 

Therefore 7 ~ 1/a for tq ^ \a\. The dimer is thus a real (virtual) bound state for 7 > 
(7 < 0), and in either case it is close to being a zero-energy resonance. 

P(po,p) has a spurious pole at high momenta (p = |p| ~ l/''"o), which is beyond the 
valid regime of SREFT and introduces unphysical short-distance effects. Here, we consider 
the physical dimer to propagate with a low momentum p ~ 7, and so expand the dimer 
propagator in powers of ro7 and r^p: 

»I'(P..P) = E -'W (^)"(o^W -mpo + pV4 )". (8) 

„ -7 + \/-^Pq + P /4 -le ^ 2 / \ J 

The n-th order dimer propagator is thus 

47r /rn\" (7 + V^-^^Po + pV4) 
mg' \2) -7 + v/-mpo + pV4 ' 
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FIG. 1: Leading-order atom-dimer amplitude, resulting from the iteration of two- and three-body 
interactions. 

Here we will consider up to n = 2 in our N^LO analysis. 

Three-body observables may be computed if we can solve for the atom-dimer scattering 
t-matrix, t. As illustrated in Fig. [T| interactions contributing to t{q,p) include the exchange 
of an atom between dimers, the three-body counterterm, and the iteration of these two 
diagrams to arbitrary order. This results in the solution to the non-relativistic Faddeev 
equation including both two- and three-body contact interactions, which is called the mod- 
ified Skorniakov-Ter-Martirosian (STM) equation [T71 [18] . The coupling constant for the 
three-body conterterm, h, is introduced to ensure that the resulting three-body observables 
are cutoff independent. 

In our perturbative analysis, we expand the quantities involved in the modified STM 
equation in powers of •yvQ and kvQ up to N^LO: 

Vipo,p) = V^'\po,p) +V^'\po,p) +V^'\po,p) + ■ ■ ■ 
t{q,p; E) = to{q,p; E) + E) + t2(g,p; E) + ■ ■ ■ 

U = Ho{A) + H,{A) + H2{E,A) + ■■■ , (10) 

where "H = K^h/2mg'^ . The on-shell t-matrix is related to the S-wave scattering amplitude 
T{k) through: 

oo 

T{k) = Zj2tn{k,k;E), (11) 

n=0 

where magnitudes of incoming momentum p and outgoing momentum q equal to the on- 
shell value k, which satisfies mE = 3/c^/4 — 7^. The renormalization factor Z can also be 
expanded in powers of •yrQ-. 

00 

Z = J2Z^, Z^ = ^^{^r,r. (12) 

n=0 ^ 
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We now review established results at LO and NLO which we use to calculate ^2, and, 
from it, three-body observables at N^LO. 



A. The Leading-Order Three-Body Amplitude 

The leading-order (LO) to is calculated via the modified STM equation, which, projected 
to the S-wave, is: 

k{q,p;E) = Miq,p;E) + - [ dq' / ' =M{q,q' ■ E)k{q' E), (13) 

Jo -7 + v3g'V4 - rnE - le 

where 

1 ( q^ + qp- mE\ 2Ho{A) 



M{q,p; E) = - log \ '\ ' — + (14) 

qp yq"^ + p"^ — qp — mE J 



The t-matrix in Eq. (13) is rescaled as to{q,p] E) = mgHo{q,p] E) to absorb the unphysical 
constant, g. 



In the bound-state region [E < — 7 /m if 7 > and < if 7 < 0), the ie in Eq. (13) 
can be eliminated, and to is therefore real. In the scattering case, to is complex. To simplify 
our calculation, we introduce a real K-matrix that satisfies a modified STM equation with 
the ie in Eq. (13) replaced by a principal- value integration (denoted by V): 

koiq,p;E) = M{q,p;E) + -P ^ dq' Miq,q'; E)koiq' ,p; E). (15) 

TT io -7 + v3g'V4 - mE 

The relation between t- and K-matrix at LO is 



ko{k,p;E) 
l-t^Ko{k,k;E) 



to{k, p- E) = .g^fc , (16) 



The fully-off-shell t-matrix at LO is related to the K-matrix (see, e.g., Ref. (SB]) by 

- f ^, , .8'fk ko(k,q;E)ko(k,p;E) 

i„fe,p;£;) = A-„(,,p;£;) + .— ^-^p-^^^. (17) 

The three-body force parameter Ho is tuned to fit one three-body observable to ensure 
a cutoff independent result. Hq has also been studied at an analytic level in Refs. [HI 135] . 
where the result 

ff^f^jy^ ^ ^ sin(go ln(A/A) + arctan(go)) ^^g^ 
sin(so ln(A/A) — arctan(so))' 
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valid up to (9(1/A) corrections, was found. In Eq. (18) sq = 1.00624, and A is the parameter 
determined by the LO renormahzation condition. For example, the relation between A and 
the atom-dimer scattering length is given by [S] 

aad= [1.46 + 2.15 cot (so In(aA) + 0.09)] a. (19) 



Ref. [18] deduced c = 1 in Eq. (18), based on an analytic calculation; while, in Ref. [35] , 
Braaten et al. found c = 0.879 yielded much better agreement (~ 10~^) with the numerical 
value of Hq. This ~ 10% effect originates from the details of the regularization. For example, 
when a hard-cutoff regularization with a finite A is performed in the numerical calculation, 
the asymptotic behavior of Kq at ~ A is distorted from the analytic result, which 



assumes A — )• oo (see Sec. VA). In order to ensure that such distortion does not affect the 
running of sub-leading three-body forces in our perturbative calculation, we first calculate 
the LO Kq at a large cutoff Aqo and insert Kq in perturbative integrations with a new cutoff 
at a lower value, A < Aqo. By doing so, the numerically calculated Kq agrees with the 
analytic expression up to g,j9 ~ A. The discrepancy induced by details of regularization is 
thus limited to 0(A/Aoo). However, the i^o needed to renormalize Kq in these calculations 



is still sensitive to the details of the regularization at Aqo and is described by Eq. (18) with 
c = 0.879, in agreement with Ref. [35] . 



B. The Next-To-Leading-Order Three-Body Amplitude 

The next-to-leading-order (NLO) calculation of effective-range corrections, ~ tq, to the 
three-body amplitude was first done perturbatively for fixed 7 by Hammer and Mehen 
in Ref. [IH], and recalculated for variable 7 by Ji et al. [^ EH J. Contributions to range 
corrections at NLO are illustrated in Fig. |2| These graphs lead to the O(ro) piece of the 
atom-dimer t-matrix, which is calculated in perturbation theory using the LO t-matrix 



obtained from Eq. (13). The NLO on-shell t-matrix, projected to the S-wave, is then 

7+ ^3g'2/4 - mE 



ti{k,k;E) 



A 



dq q - , 

-7 + v/3g'V4 - niE 



le 



--pQ{k,q';E) 



2^1 (A) 



A2 



2 o'2 

1 + - / dq' , ^ 

TT Jo -7 + v3g'V4 - mE - ie 



io{k,q';E) 



(20) 



where ti and Hi are rescaled as ti = VQing'^ti, and Hi = tqHi. Here we only show the elastic 
scattering case p = q = k, which is what we are interested for the computation of, e.g. the 
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FIG. 2: The diagrams for the atom-dimer t-matrix at NLO. The NLO piece of the dimer propagator 
is denoted by the square labeled "1" and NLO corrections to the atom-dimer contact interaction 
are indicated by the circles labeled "1". 

phase shift. 

Using a similar approach to the LO case, the complex t-matrix at NLO, ti, can be related 
to a real K-matrix, which satisfies an expression with principal- value integrations [23]: 



g2(7+ ^3g74 - mE ' 



K,{k,k;E) = -V / dq 

7^ Jo -7 + V3gV4 - mE 



kl{k,q-E) 



2^1 (A) 
A2 



1 + -V 



TT 



dq 



q^Ko{k,q;E) 
-7+ y/3q^/A-mE 



(21) 



The NLO pieces of the on-shell t-matrix and K-matrix are related by 



ii{k,k]E) 



1 - i 



Ko{k,k;E) 



k,{k,k;E) + i^ko{k,k;E) 



(22) 



The introduction of a real Ki not only simplifies the calculation but also straightforwardly 
leads to the NLO effective-range correction to the phase shift. 

Refs. [221 [23] pointed out that the three-body force Hi{A) will in general be 7-dependent: 



#i(A)=A/iio(A)+7/^ii(A), 



(23) 



where both hiQ and hu are functions of A, whose analytic expressions were derived in 
Ref. hiQ is of 0{A^) and is log-periodic in A 

37r(l + si) a/1 + 4sg - cos (2so ln(A/A) - arctan 2so) 



/iio(A) 



(24) 



64-^1 + 4so sin^ (sq ln(A/A) - arctan sq) 

The 7/2.11 piece of Hi is ~ 7 In A (see Ref. [23] for details). If we are only interested in 
the fixed-7 case, as is relevant for three-nucleon systems or the case of ^He trimers, then 



it is appropriate to think of Hi as one overall number, i.e. the 7/iii part of Eq. (23) can 
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(A) 

FIG. 3: N^LO t-matrix dia gram part A: Insertion of N^LO dimer propagator (square labeled "2"). 

be combined with hiQ. The two pieces can, however, be disentangled by experiments in 
systems with a variable scattering length, such as the recombination of ultracold ^Li or 
Cs atoms [371 EH]. 



III. THE THREE-BODY SCATTERING AMPLITUDE AT N^LO 

We now calculate the next-to-next-to-leading-order (N^LO) effective-range corrections 
to the atom-dimer amplitude perturbatively. We focus on the elastic-scattering situation 
p = q = k and study N^LO corrections to three-body bound-state energies and the S-wave 
elastic-scattering phase shift in three-body systems with fixed scattering length (i.e. fixed 
7) and effective range. 

Since three-body energy is conserved in an elastic channel, we omit the third argument, 
E, in expressions of t-matrices in this section to avoid cluttering notation (i.e. to{q,p) = 
to{q,p;E)). We will restore the third argument in subsequent sections. 



A. Diagrams and t-matrix 

In our perturbative analysis at N^LO, we insert effective-range corrections to the dimer 
propagator and the three-body counterterm between LO t-matrices, which are obtained from 



the modified STM equation (13). Compared to the calculation of the t-matrix at NLO, here 
we insert both NLO and N^LO terms and consider their contributions to the 0{r'^) piece 
of the atom-dimer t-matrix. The resulting diagrams are categorized in classes from A to E, 
that are illustrated in Figs. [3jj8| 

Fig. shows the effect of the N^LO correction to the dimer propagator. This constitutes 
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(B.l) 



(B.2) 



(B.3) 




(B.4) 



FIG. 4: N^LO t-matrix diagram part B: Insertion of N^LO 3-body contact interaction (circle 
labeled "2"). 

diagram A, whose contribution to the N^LO t-matrix is 



<(k,p) 



V^^\qo, q) iS{E - go, -q) ito(k, q) zto(q, p) 



(25) 



where we only consider the on-shell case that |k| = |p| = k and mE = 3A;^/4 — 7^. By 
projecting t^(k, p) onto the S-wave we obtain 



t^(A;, k) 



2Tcmg'^ 



— 7 + y 3g2/4 — rni^ — ze 



tl{k,q). 



(26) 



At N^LO three-body (= atom-dimer) contact interactions will be required to renormalize 
the amplitude. We denote these as H2. Based on the power-counting arguments given in 
Ref. pE] we anticipate that, in addition to being dependent on A, H2 will also depend on 
the three-body energy E, i.e. we have H2{E,A). Fig. |4] defines class B, the contribution of 
such contact interactions to the atom-dimer scattering amplitude. From the Feynman rules, 
it's formulated as 



.2mg^H2{E,A) 
' A^ 



1 + 



dS 



iV^'\qo, q) tSiE - go, -q) ^to(k, q) 



/^4 
iV^"\qo, q) iS{E - go, -q) «to(p, q) 



Here we define the part in the square bracket in Eq. (27) as a function £(p): 

£(p) = 1 + r 7^ tV^'\qo, q) tSiE - go, -q) 2to(p, q), 
Jo (27r)* 



whose S-wave projection is 
Cik) = 1 + 



dq 



-7 -|- ^J'iq^ 1^ — mE — ie 
13 



to{k,q). 



(27) 



(28) 



(29) 




(C) 

FIG. 5: N^LO t-matrix dia gram part C: Double insertion of NLO dimer progators (squares labeled 
"1"). 

The function L will also appear as part of the N^LO contributions from diagrams in classes 
D and E. The S-wave projection of in the on-shell case is then 

fMAKU) = ^'"^^gf.^) ^^,,). (30) 

This accounts for contributions directly from the N^LO dimer propagator and N^LO 
three-body counterterms. We now compute contributions arising from two insertions of 
NLO terms. These constitute classes C-E. 

Fig. |5] represents double insertions of the NLO part of the dimer propagator: class C. 
The resulting contribution to the N^LO t-matrix is 

xito(k, q) ito(q, q') ^^o(q', p), (31) 
whose on-shell S-wave projection is 



A A 

,2 



^2^2^4 J J ^3gV4 -mE -it -7 + v^3g'V4 - mE 



le 







Xto{k,q)to{q,q')to{k,q'). (32) 

Fig. [6] represents class D: graphs with two insertions of the NLO three-body counterterm. 
The on-shell contribution can be re-expressed in a form that factorizes counterterms and 
closed loops (see Fig. [7]for a diagrammatic formulation). The loops Xi and X2 defined in 
Fig. [7] are 

^Il = l^.j (^^^^^°Hgo,q)^5(i5-go,-q)^to(q,q') 

xtV^'\q'„c^)tSiE-q'„-c^) , (33a) 
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(B.l 



fD.51 








fD.6l 




fD.7l 




fD.81 



FIG. 6: N^LO t-matrix diagram part D: double insertion of NLO 3-body contact interactions 
(circles labeled "1"). 



X 



[ (DY X [ 1 + 





^0 




FIG. 7: An alternative way of writing part D of the N^LO t-matrix diagrammatically. 

= I ^^zV('\q,,q)zS{E-qo,-q) , (33b) 
which are both only functions of the three-body energy. In terms of these the on-shell S-wave 
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^0 



(E.l 




E.2 




(E.3) 




E.4) 



^0 




E.5 




E.6 





FIG. 8: N^LO t-matrix diagram part E: These arise from the insertion of one NLO dimer propagator 
(squares labeled "1") and one NLO three-body contact interaction (circles labeled "1"). 

projection of class D is 



^2 (^5 ^) 



2m^2^i(A) 



C\k) [X1+X2] 

-A 



7rA4 

A A 



-7 + -\/3g^/4 — mE — ie 



+ 



Timg- 



J J dqdq' — 



J2 







7 + V3g74 - mE — ie —7 + 



(34) 



Fig. |8] includes terms with one insertion of the NLO dimer propagator and one of the 
NLO three-body counterterm. These make up class E. Their factorization is shown dia- 
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2[©] X [1+ (t^l X [^^ + ^^^1 

FIG. 9: An alternative way of writing part E of the N^LO t-matrix diagrammatically. 
grammatically in Fig. |9| where the loop integrations J3 and 1^ are, respectively, defined as 



f (Tq 

= iWfl 0ji''^^'^M^S{E-qo,-ci)zV('\qic^)zS{E-qi-ci') 

xzio(q,q')^io(q,k). (35b) 

Therefore, the contribution of class E to the S-wave on-shell N^LO atom-dimer t-matrix 

is 



^2 {k, k) = — C{k) [X3 + X4J 



7rA2 







dq /q 2 /a j^ =^o{k, q) 

—7 + a/ oq^/A — niE — %e 



A A 



H / / ctgdg , , to[k,q )to[q,q ) 

nrng^ J J -7 + V3gV4 - - «e -7 + 73g'V4 - - 







(36) 



B. The atom-dimer K-matrix at N^LO 



The integrals in the previous subsection include terms ie in the denominator, and so, in 
the scattering case, the N^LO t-matrix t2 is complex. Each term in the dimer propagator's 
Tq expansion ([o]), has the same denominator, and so they all have a pole at the on-shell 
momentum for atom-dimer scattering. We therefore separate the propagator into a Cauchy 
principal- value and imaginary part: 



-7 



y^3qyI^~mE^^ -7 + y/3q^/A-mE 3A; 
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In the following exposition, we do not include the principal-value symbol. Except where 
especially noted, all the integrals with a real-number singular pole should be considered as 
principal-value integrals. 



We now isolate the delta-function part in Eq. (37) for each piece of the N LO contribution 



to the t-matrix. In a bound-state problem with a negative three-body energy E, these 
imaginary parts need to be omitted. It is therefore useful to relate t2 to an N^LO K-matrix, 
in which only real integrals appear. 



After using Eq. (37) in the expressions for t^-tf, and employing Eqs. (17) and (21), we 
find that the full expression for the N^LO piece of the t-matrix takes the form 

K2{k,k) , ^g^},{jKo{k,k) + K,{k,k)' 



t2ik,k) 



l-i^Ko{k,k) 



l-i^Ko{k,k) 



(38) 



where we have rescaled t2 according to t2 = r^mg'^t2- Here 



K2{k,k) 



1 

2^ 







-7+ ^Sq^/A-mE 



2H2{E,K) 



A A 



A 



dq 



q^Ko{k,q) 



+ - 




dqdq 



7rA4 







TT 



7+ ^3gV4-mE 

, + y/Sq^/A - mE ) q'^{-f + ^JSq'^/A - mE ] 
-7 + ^/3q^/A - mE -7 + ^/Sq'"^ /A - mE 



xkoik,q)koik,q')ko{q,q') 
q^ko{k,q) 



-7+ y/3gV4 - mE 



dq 



-7+ ^3q^/A-mE 



1 + 



dq' 



q"Koiq,q') 



-7+ y/3q'^/A-mE 



4^1 (A) 
7rA2 



TT 



q'ko{k,q) 



-7+ y^3(f/A^^mE 



dq 



q^{n+ ^/'iq^lA - mE 



, -Ko{k,q) 1 + 

-7+ V3g2/4-mE ''^'^M tt 



dq' 



q"Koiq,q') 



-7+ ^3q'^/A-mE 



+ K!^ + Kli + K!f + K. 



2 ) 



(39) 



with the N^LO three-body counterterm parameter rescaled as H2 = '"o-^s- Each of the 



integrations summed up in Eq. (39) is identified with K2, n = A-E, thus associating 
each contribution to the N^LO K-matrix with a particular category of diagrams defined 
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in Subsection III A Since Eq. (38) can also be derived by considering the general relation 
between t- and K-matrices and expanding to second order in tq the fact that we obtain 
t2 with this structure is a non-trivial check that we have perturbative unitarity, i.e. our 
t-matrix satisfies unitarity up to corrections of 0{rQ). 



IV. THREE-BODY OBSERVABLES AT N^LO 

Here we discuss how to calculate three-body observables at N^LO within our perturbative 
approach. In an atom-dimer scattering state, we consider the elastic-scattering phase shift 
in the S-wave; for trimer bound states, we focus on the binding energy. These quantities are 
expanded in a series of and/or krQ, with the leading order representing universal physics. 
We derive NLO and N^LO shifts in these quantities as corrections beyond universality. 



A. S-Wave Phase Shifts 



The S-wave amplitude T{k) of an atom-dimer scattering state is related to the phase shift 



by 



T{k) 



Sir 



(40) 



m k cot 6{k) — ik 
k cot 6 can be expanded in powers of tq 

k cot 6 = k cot 6o + ro[k cot S]^ + rl[k cot 6]^ + ■ ■ ■ , (41) 

where []-^ and [jg indicate respectively the NLO and N^LO part of kcot6. Therefore, the 
scattering amplitude T{k) is expanded correspondingly as 

[k cot 6] ^ 



T{k) 



37r 
m 



fccot^o-^^ {k cot 5o- iky 



{[k cot 5] J 



[k cot 6]r 



{k cot 6o — ik)"^ {k cot 5q — iky 
= To{k)+T,{k)+T2{k) + -- - . 



+ ■■■ 



(42) 



The LO amplitude is related to the K-matrix by Eqs. (11), (12), and (17): 

Toik) = Zoto{k,k;E) 
^''\{k,k;E) 



m 
3n 



ik 



(43) 
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which means that we recover the LO relation between phase shift and K-matrix 

3 ~ 

fccot^o = —KQ\k,k;E). 
87 

Similarly, the NLO amplitude is related to K-matrices as: 

Ti{k) = Zito{k,k;E) + Zotiik,k;E) 

= ro— [-fio{k, k- E) + h{k, k; E)] , 
m 

which leads to the relation for the NLO part of k cot 6: 

[kcot6]^ = -^k^\k,k-E) (ki{k,k-E)+-fko{k,k;E) 
87 ^ 



(44) 



(45) 



(46) 



Proceeding in the same way at N^LO, we derive the N^LO relation between the scattering 
amplitude and K-matrices as: 



T2{k) = Z2toik,k;E) + ZMk,k;E) + Zot2ik,k;E) 



, 37r 



m 



87' 



+ 



87' 



K,^K2-jK^-Kf/Ko 



87-^0 



(47) 



where we substitute the expressions of ti (22) and t2 (38) into Eq. (47). After recognizing 



that part of the expression (47) can be written as /ccot^o and [A; cot 5] 1 we compare it with 
the N^LO piece of Eq. (|42|, and find 



3 ~ 



[k cot 5]^ = K^\k, k- E) K2{k, k; E) - -fKi{k, k; E) - ^ 



kf{k,k;E) 



Koik,k;E) 



(48) 



B. Three-Body Bound States 

When the three-body energy is below the atom-dimer threshold {E < for 7 < and 
E < — 7^/m when 7 > 0) all integrals become purely real and the principal- value prescription 
can be dropped, since no singularity appears in the integrand. Thus i = K in this case. 

A three-body bound state exists with the binding energy B when the three-body K- 
matrix has a pole at E = —B. The K-matrix can be expanded around the position of its 
singularity 

(49) 



Kiq,p;E) = ^^+niq,p;E), 



E + B 
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where functions Z and 7^ are defined as the residue and regular part of this expansion. 
In our perturbative approach, the binding energy B is expanded in powers of tq as 

5 = Bo + roBi + r2S2 . (50) 



We can also expand each term in Eq. (49) in powers of tq: 



E + B^ + roB^+rlB^ 



+ 7^o + ro7^l + rln^ 



(51) 



The next-to-leading-order energy shift is then 

Bl = 



(52) 



which is independent of incoming and outgoing momenta, q and p, so we can take q = p = k 



and calculate Ki{k^ k; E) from Eq. (21) just for convenience. In fact, the residue function 
Z(){q,p) takes a separable form [23j: 



Zoiq,p) = T{q)T{p) 
with r(g) satisfying a homogeneous integral equation 

r(g) = I [\q'Miq,q';-B,) 

^ Jo -7 + V 3g'V4 + mBo 

Therefore, Bi can, instead, be calculated from 



(53) 



12 



r(g')- 



(54) 



5i 



1 



7r 



-7 + V3gV4 + mBo 



(vrA)^ 



A 



g'r(g) 



-7 + V3?V4 + rri^o 



(55) 



Similarly, the N^LO shift of the binding energy is related to residues of the double pole 



and triple pole in the expansion of Eq. (49): 
Bo = 



{E + BofK^{q,p- E) + {E + Bo)B,K,{q,p- E) 



Zo{q,p) 



(56) 



In Eq. (56) {E + i?o)^-f^2 has a pole a.t E = —Bq which is canceled by the corresponding 
pole of {E + Bq)BiKi. This cancellation is seen explicitly if we derive an expression for B2 



analogous to Eq. (55) for Bi. 
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First, we insert Eq. (53) into Eqs. (21), and expand {E + Bq)Ki as 

Jo -7 + A/3aV4 + m5n 



r\k) (A) 



dq 



-7 



-7 + \/3gV4 + m^o 



+ v/3g74 + m5o 



g^7^o(fc,g;-i?o) 

rfg , ^= 

-7 + V3gV4 + mBo 



2r(fc) 



TT 



/•A 



2 7+V3gV4 + m5o 

, , , =^(g)7^o(fc,g;-^o) +•••, (57) 

-7 + V 3^' /4 + mBo 

), that vanish in the hmit E — )■ —Bq, and the 



where elhpses indicate terms of O {{E + Bq)^ 
pole at E = —Bq is now exphcit. 

{E + BqYK2 can be computed using Eqs. (21), (39), and (53) in a similar way, and 
expanded around E = —Bq. It too has a first-order pole at = —Bq. Inserting {E + 
BqYK2S and {E + Bq)Kis pole expansions into Eq. (56) and replacing Bi by the expression 
(55), we find that the terms singular ai E = —Bq cancel, leaving: 



B2 = / dqq^ 

8iJ2(-5o,A) 



t±VMI±J!}?^Y\q) 
tuBq 

q^T{q) 



2 — / v , = 

2vr Iq -7 + 73gV4 + 



(ttA)^ 



-7 + \/3gV4 + mBQ_ 



A A 



~Z:2 / / dqdq q q ====== ======1 {q}L {q } l<,Q{q, q , -Bq) 

^ V3gV4 + uiBq -7 + V3g'V4 + tuBq 





7r3A4 



-7 + 



8^1 (A) 



(ttA)^ 



Jo 

r 



^A 



v/3gV4 + mBo 



-7 



-7 + v/3g2/4 + m5o I ^ ^ tt 



, q''nQ{q,q';-BQ) 
^ -7+ V3g'74 + m5o^ 



g- 



'r(g) 



(ig , 

-7 + V3g2/4 + m5o 



2 W^^^r(g) (l+'-rdq'^t^M=^ 
Jo -7 + V3gV4 + uiBq ' ' V vr io -7 + \/3g'74 + ^ 



-7+ V3g'74 + m5o^ 



(58) 



which is finite at £" = —Bq. This form for i?2 shows that, after renormalization in which 
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any significant cutoff dependence in B2 will be cancelled, B2 is only a function of the LO 
binding energy Bq. It is independent of incoming and outgoing momenta. 



V. ASYMPTOTIC BEHAVIOR OF LO K-MATRIX 

We list several important asymptotic features of the LO K-matrix, Kq, in this section, 
for both half- and fully-off-shell cases. These properties are essential in the calculation 
of H2{E,A) when the N^LO renormalization is performed. A detailed derivation of these 
asymptotics is given in Appendix \K\ 

General features of Kq^s asymptotic behavior are renormalization-condition independent. 
However, the normalization factors of Kq, in both half- and fully-off-shell deter- 
mined by a specific renormalization condition at leading order. To compare analytic forms 
of -ft'o's asymptotics with numerical values, we employed a specific physical condition in this 
section: the atom-dimer scattering length is fixed to 'yaad = 1.5 at leading order. Fitting the 
normalization factors and a few other parameters to their numerical values in this physical 
choice, we can study general features of Kq^s asymptotics with high accuracy. 



A. Asymptotics of half-shell K-matrix KQ{k,p; E) 

The asymptotic behavior of the half-on-shell Ko{k,p;E) has an expansion in powers of 
7/j9 and mE/p^ at large p: 

koik,p; E) oc (1 + + + + ■ ■ ■ ) ' (59) 

where ellipses indicate higher-order terms of 0{'y^/p^) or 0{m?E'^ /p^). The constants 
and Ci are computed to be (see Appendix of Ref. [20] where we have corrected their result 
for Di by a factor of 1/3) 

„ ( 2Y I(isn-n) 

C. ^ , (60b) 

1 - I[isq - 2) 
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where I{s) and L{s) are functions calculated from Mellin transforms (see Ref. |59] for deriva- 
tion): 



lis) 



Lis) 



8 sinf 
VSs cos ' 



8 sin 



2 

7r(s-l) 



COS 



(61a) 



(61b) 



After fixing a renormalization condition at LO, the K-matrix is a real number associated 
with observables. Therefore, we can rewrite its expression as 



with 



Ko{k,p; E) = ay{k) 
<i>n{p) 



1 ^ ^ TYlE 

-Mp) + -^Mp) + -^Mp) + ^^i(p) + 
p p^ p-^ p-^ 



\Dn\ sin ^so In ^ + arg 

/ p 
I Ci I sin l^so In + arg C 



(62) 

(63a) 
(63b) 



The normalization factor a^{k) in Eq. (62) is generally a function of both 7 and k. Since in 
this paper we only consider a system with a fixed scattering length, 7 is also fixed and so is 
treated as a constant parameter. With this constraint in our analysis, only varies with 
the on-shell momentum k. 

Furthermore, at the zero-energy threshold, = 0, the asymptotics of KQ{k,p; E) can be 
expanded just in powers of 7/p, as 



Ko{k,p;E) = a^{k) ^ 



n=0 



7 I I f P 

Dn\sm (^soln^ + argL'„ 



(64) 



where k = ^ at E = 0. 



B. Asymptotics of fuUy-off-shell K-matrix KQ{q,p;E) 

When the incoming and outgoing momenta are both large, i.e. q,p^ k,'j, the fully-off- 
shell K-matrix KQ{q,p; E)^s asymptotic behavior is dominated by the leading piece, ~ l/{qp), 
in the expansion at large p and q. Here, we focus on the behavior of this leading term and 
do not consider sub- leading terms in the expansion. 

Using power-counting analysis and numerical fitting (see also Ref. [3^) we find that 
Ko{q,p; E) at large momenta p, g ^ /c, 7 is dominated by 

ko{q,p;E) = ^[J^,(q,p) + g,(q^p;E)] . (65) 
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The energy- independent part of Kq is represented by the function J'q, which obeys 

Mq,p) = { _ , (66) 

b-fM(l)Mp) + PiP^Q)^ if (p > g) > fc, 7 

where is a constant since 7 is fixed (its value is determined by the LO renormahzation 



condition), 0o is defined in Eq. (63), and 0o is defined as 



Mp) = cos (^solnjj . (67) 
The energy-dependent part of Kq is written as a function Qq: 

go{q,p;E) = c^{E)MQ)Mp) , (68) 

where the factor c^(-E') only depends on E at fixed 7, and is also determined by the LO 
renormahzation condition. 

The product 0o0o iii function J-'q, together with 0q in function Qq, dominate off-diagonal 
(i.e. p ^ q OT p q) elements of the LO K-matrix: analysis of the modified STM equation 



( 13 ) shows that p tends to zero if q ^ p or p ^ q. The function p thus represents a remainder, 
which plays an important role in the near-diagonal (i.e. p ^ q) K-matrix elements. The 
remainder function depends on the renormahzation condition and we have only determined 
how to numerically fit it to a chosen parametrization. 



Here we will verify each piece of the approximate form of the fuUy-off-shell Kq in Eqs. (65 



68) by comparing these results with numerical calculations. In these calculations, we first 
fix the LO atom-dimer scattering length to be 70^^ = 1-5, and then calculate Kq at = 0. 
Firstly, at = 0, we can combine functions Tq and Qq and express Kq as 

~, U7 0o(g)xo(p) + p(g,p), if(g>p)>7 

^o(g,p;0) = < , (69) 

U7Xo(g)0o(p) + p(p,g), if(p>g)>7 



where xo obeys 



Xo{p) = + c^(0) sin {^Q In ^ + arctan ^r^) • C^^) 



We show in Fig. 10 that the first term, proportional to </>o(q')xo(p) or Xo(g)0o(p)5 in our 



approximation (69) agrees well with the numerical results when q and p are not close to 



each other. Here, in order to compare with Eq. (69) for Kq, we performed a best fit to 
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10^ 10^ 10^ 10^ 10^ 10^ 10^ 

p[y] 



FIG. 10: Amplitudes qpKo{q, p; 0) as functions of the momentum p at fixed values of q: q = 2x 10^7 
(upper panel, a) and q = 2 x 10^7 (lower panel, b). Dots are the numerical results. Solid lines 



(red) represent the analytic function defined by the first term of Eq. (69). In each case, Hq{A) is 
chosen to reproduce ^Uad = 1-5. 

the amplitude and phase of the function xo within an accuracy of 10^'^, which determines 



= —2.012 and c^(0) = —8.612. Meanwhile, A is parametrized in Eq. (19), and, when 
fitted, turns out to be A = 4.4217 O.I57 correction to the result from Ref. |10], to obtain 
better accuracy). 

The remainder function p{q,p) describes discrepancies between the numerical values of 
KQ{q,p;0) and this off-diagonal part (i.e. 4>o{q)xo{p) or Xo{q)(Po{p))- Since the final ex- 
pression for KQ{q,p;0) displays discrete-scale invariance we choose a form for p{q,p) which 



depends only on the ratio of q/p. And indeed, we show in Fig. 11 that p{q,p) depends only 



explicitly on the ratio q/p if q/p ~ 1 (the same parameters used in Fig. [lO| were chosen in 



producing Fig. 11). We must also respect the symmetry under the exchange of incoming 
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FIG. 11: Remainder p{q,p) as a function of the ratio q/p, with each curve representing a different 
value of p. As in Fig. [lOl the results are for a renormahzation condition 70^^ = 1.5. 



and outgoing momenta {KQ{q,p; E) = Ko{p, q; E)). Hence we fit p to the approximate form: 

q/p, if q < p 

(71) 

p/q, ii q > p 

where a and (3 are renormalization-scheme dependent. The particular function employed in 



p(x) 



axe 



X 



Eq. (71) has no theoretical justification and is chosen solely for its ability to give a good 
fit. Once the same parameters as in Fig. 10 are chosen, we find a = —1.665 and /3 = 1.322, 



which are fitted at the region x ~ 1 within an accuracy of 10" 



In Fig. 12, the numerical result of p{q,p) is plotted as a function of q/p at a fixed value 
of p, and it agrees well with our analytic formula with the parameters given above. As 



suggested in Fig. 11, such an approximation loses accuracy if q/p is not close to 1. We 
show in Sec. VI B that the effect of this on the determination of the N^LO three-body-force 



parameter, h2o, is small. 



Furthermore, we can also obtain c^(£')'s value with simple algebra based on Eq. (65): 

ko{q,p; E)-ko{q,p;0) 



c^{E) = c^(0) + 



0o(g)0o(p) 



(72) 
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FIG. 12: Remainder p{q,p) as a function of the ratio q/p, where p is fixed to 10^7. The numerical 
result (dots) is renormalized to reproduce 'jUad = 1-5. The analytic expression (red solid line) is a 



fit (71) to the numerical result. 



Assuming our approximation is correct, we should observe that is q- or p-independent 



when we insert numerical results for Kq into Eq. (72). The outcome of this test is shown in 



Fig. 13, where we plotted the second term of Eq. (72) at different values of q, p and E. It 



shows that Eq. (72) is a function of E alone: it does not change when g or p varies, provided 
that q > 10^7. The divergence below this value is due to higher-order corrections to Kq (i.e. 



0(7/^,7/]?)), which shift the numerator in Eq. (72) from zero at (j)o{q) = or 0o(p) = 0. 

We also expect that the near-diagonal part of Kq, p{q,p), is energy independent when 
only the dominant part of Kq is considered. Our expectation is supported by our findings 
in Fig. 13, since no bump is observed in the curves at q = p = 10^7 or q = p = lO'^'y. 



In this way, we describe the asymptotic behavior of the dominant part (~ l/{qp)) of the 
fully-off-shell Ko{q,p; E). 
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FIG. 13: c^{E) — c-y(O) as a function of q at different fixed values of p and E. When p is fixed to 
10^7, three lines respectively represent the cases for mE = —7^/3 (red solid line), mE = —27^/3 
(blue dashed line) and mE = —7^ (orange dotted line). When p is fixed to 10^7, dots represent 
the cases for mE = —7^/3 (circles), mE = —27^/3 (squares) and mE = —7^ (diamonds). Results 
are renormalized to reproduce 70^^ = 1.5. 

VI. CUTOFF DEPENDENCE AND RENORMALIZATION 



In the absence of additional renormalization, K2 has cutoff dependence ~ at N^LO, 
which renders unrenormalized results from N^LO calculations unphysical. A sub-leading 
short-range piece of the three-body/atom-dimer interaction, here denoted H2, is thus needed 
at N^LO to absorb this cutoff dependence. Since H2 contains powers up to (9(A^), we write 
the N^LO three-body force parameter H2 as 



H2 = A^h2o{A) + mEh22iA), 



(73) 



In this section, we study both the leading piece of if 2, ^20, and its sub leading piece h22- 
7-dependent pieces of this three-body force are absorbed in /120 as a fine-tuning modification 
in problems with a fixed-but-finite 7. Renormalization must be carried out with such terms 
made explicit if we want to study systems with varying 7 
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A. Relevant integrals 



The principal- value part of the integration in C{k), with replaced by Kq, appears 
repeatedly in the LO, NLO and N^LO perturbative calculation, i.e. in Eqs. (15, 21, 39). 



Here we denote the integral itself (with the infrared regularization dropped) by ^: 



7r 



(74) 



-7+ ^3q^/4:-mE' 

where, as usual, a principal-value prescription is to be understood as the means of dealing 
with any pole that is present in the integrand. Upon inserting the asymptotic form of 



Ko{k,q), (64), up to terms ~ l/q^, we then have: 



27 472 2mE 

dq{l + — + — + 

y/'iq 3q^ 3 q^ 



^0 + -01 + ^ 
q q^ 



^ mE 

2 + ^r^i 



ysvr ' ' J "V V3 ? -iq" ^ y V q q^ q 

In order to simplify the notation, we define integrations that contain functions 



.(75) 



as 



We can also define 



!^m(A) 



1 

vr 
1 

TT 



dq q 



m—l , 



dqq^'-^tjji. 



(76) 



(77) 



Since 0„ and ipn are C(l), the subscript m of the defined functions ^m"* and indicates 
they have cutoff dependence, ~ A™, except that m = indicates the possible presence of 



both terms ~ InA and 0{1). Eq. (75) can then be written as an expansion in 7/A and 
mE/A^ 



4:a-y{k) 



(1) 



3 



(2) 



+mE 



s(0) 



(78) 



The infrared regularization of some integrals (e.g. (pf'\ ^q"'* and ^0^'') in (78) could affect 
the 0{1) pieces of C. However, we have numerically verified that, when computing C{k), 



the combination of infrared parts of these integrals is canceled by the "1" term in (29). This 
cancelation yields C{k) = ^. 



When 7 is fixed, the 7-dependent terms in Eq. (78) can be excluded in the analysis of 



cutoff dependence, since they can be absorbed by fine tuning of Hq, hio and /i2o respectively 
in the LO, NLO and N^LO renormalization. In the rest of this section, we do not consider 
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such 7-dependent terms. Therefore, of all the functions 0„, only 0o is needed for our purposes 
here. We can thus further simplify our notation for (Pm^ by dropping the superscript (n): 



1 '•^ 



<Pm{A) = - / dqq'^''(j)o. (79) 



^'s expression is then reduced to 

(k■^ r /9 M 

(80) 



^ _ 4a^(fc) 



where each integral in Eq. (80) is calculated in Appendix [B| <Pi is used for calculating h2o, 
and the addition of terms ~ mE is needed for the calculation of h22- 

For later convenience, here we also define some other integrals that are needed in the rest 



of this section. Integrals containing 0o are defined similarly to $„: 

^„(A) = - t dqq-~^Uq) ■ (81) 



Analytic expressions for the ^„ and ^„ are given in Appendix |b] 

Integrals that contain a product of two functions, e.g. 0o and 0o; are defined by 

W^0,n(A) = I f dqq--'Uq)Uq). (82) 

where the subscript n also indicates Ws power-law dependence on A. Analytic expressions 
for U^^^n and W^^^n are calculated in Appendix |Bj 

The double integrals are defined in a similar way, for example: 

I 

y^4>m,m+n{^) = - J dqq'^-^4>l{q)$n{q) 

= ^ f dqq^-'ct^liq) J' dpp^-'Mp), (83a) 
>V<^(00),r«+n(A) = - / dqq"" ^(t)o{q)U^^,M 



1 



^ dqq^-'Uq) r dpp''-'Mp)Mp). (83b) 



where m + n denotes the integrals' overall powers of A. In Appendix [B| we write the 
expressions for and >V<^^(^^)^i+i, which are all shown to be 

~ A2. 

Another type of double integral, that contains the near-diagonal function p{q,q'), is de- 
fined by 

A A 

n^p^{A) = ^J J dqdq'cl){q)p{q,q')cl){q'), (84) 
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where the g-dependent (g'-dependent) function appears on the left (right) side of p in the 
subscript. These integrals will be calculated directly in the next subsection. We also will 
use the notation Vtp^, Q^p and Qp. If does not appear in the designated location in a. fl 



subscript then it is to be replaced by "1" in the integrand in Eq. (84). 



B. The dominant cutoff dependence and h 



20 



We now calculate the dominant cutoff dependence in each of the terms in Eq. (39) cate 



gorized from A to E. Inserting the leading piece of Eq. (62), we see that the first term, (A), 
has dominant cutoff dependence: 



= aiik) — I da ad'i(Q) = aiik) —U^^,(A). (85) 
After inserting the dominant pieces of ^ from Eq. ([80|) and the dominant piece of H2iE, A) 



form Eq. (73) we find that the second term in K2, K2 (the part of K2 that results from the 



insertion of the N^LO 3-body force) behaves as 

,2 



K. 



^-^^^^e^ = <(fc)fMA)^?(A). 



(86) 



By inserting the off-shell KqS asymptotic form (|65j), we find that the third term from 
two insertions of the NLO piece of the dimer propagator (C) is calculated as 

= ^ / / dqdq'Mq)Mq')[J'o{q,q') + Go{q,q')] 




dq4>l{q) 



+ 



aUk) 



A A 




dqdq'(t)Q{q)(t)Q{q')p{q,q') 



A rq rA /-A 



o(g) / dq (t)Q{q)(pQ{q) + / dq(f)o{q)M(l) I H (fliq" 



= a'^ik) [c,{E)Ul^^,{A) + 2b,W^^^^^)^,^,{A) + n^p^{A)]. 
where we used the fact that dq dq' = dq' dq in the last step. 



(87) 
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The fourth term (D), from double insertions of the NLO piece of the 3-body force, is 



K. 



D 



32MA) 
3A 



q + 



4 



y/SlT J V V^TT 



dq' [Mq.q') + Qo{q,q')] 



2b. 



a5(A;)c??(A)|^| qdq + c,iE) dq<Pl{q) 
dqMq) j dq'(j)o{q') + ^ j j dqdqp{q,q') 



32MA) 
3A 



c,{E)<Pl{A) 



in 



+ 26,>V^(^),i+i(A) + r],(A) 



(88) 



where the leading piece of Hi, A/iio(A), is inserted in the calculation. 

Finally, the last term (E), from the insertion of one NLO dimer and one NLO 3-body 
force, is 



^2 



4A/iio(A) a^{k) 



A2 



TT 



dqMq) {q + 



7T 



dq' [Mq.q') + Goiq,q')] 



3A 



4:71 



^^^''^^^-a'^{k)<Pi{A){^ I dqqMq)+ 



c,{E) 



dqcpliq) 



'Mq) 



+ - 



dqMq) / dq'p{q,q') 



+ - 



64/iio(A) 



dqMq) I dq'Mq) + I dqMq) J dq'Mq')Mq') 

a'Jk) i(A) <; c,{E) ^i(A)W^^,i(A) + ^^^(A) 



3A -7V-/ — V--/ I -'V-/ — v--/--V's-.-v--/ ■ 4 

+^'7["^</.</.(0),i+i(A) + >V<^{<^^),i+i(A)] + n^p{A) } . (89) 

The infrared regularization of each integral in K2-K2 can in fact contribute to the 0{1) 
part of K2; however, their contributions are two orders smaller in powers of A compared to 
the dominant cutoff-dependent pieces of K2. We therefore drop the infrared regularization 
in each integral for our convenience. By doing so, we treat expressions for these indefinite 
integrals (calculated in Appendix |B]) as exact results, with unsolved infrared pieces absorbed 
in the renormalized K2. 

In order to simplify the calculation, we regroup integrations 



+ + + = a'Jk) (Se + + ^2D + Sp), 



(90) 
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where T,e indicates the summation of terms proportional to c^{E), T,io sums up terms that 
contain a ID integral, T,2d denotes the summation of 2D integrals W, and Sp refers to the 
summation of integrals including p(g, q'). 

We first sum up the three terms oc c^{E) (in K2 , K2 and K2) by 



T.E = c^{E) W^^,i(A) + 



3A ' 



(91) 



The A-dependence here is canceled in the NLO renormalization, because hi^ is tuned at 
NLO to make the ~ A pieces in the square bracket cancel with each other: 



3A 



^?(A) 

TT (1 + Sg) + 4so - cos (2so ln(A/A) - arctan 2so) 
2A ^1 + sin^ (soln(A/A) - arctan sq) 



(92) 



which agrees with /;,io(A)'s expression (24) that is obtained in 



By substituting Eq. (92) for /iio, we write the summation of ID integrals (in K2, 
and K2) as 



4 



(93) 



^i(A) ' <Pl{A) 271 _ 

Since the subscript of all integrals defined above displays the order of their cutoff dependence, 

readers can easily verify that each term in ~ A^. 

We then sum up the three 2D integrations of the off-diagonal amplitude (in C-E), and 

obtain 

S2D 



2K 



(A) + ^^^i^w 



^i(A) 



1+1 



(A) + 



1+1 (A)]. 

(94) 



Analytic expressions for the four W's in Eq. (94) are given in Appendix [Bj Each term in 
S2D ~ A^. 

The summation of integrals containing the near-diagonal remainder function p results in: 



= n,,,iA) - 2%i^(^,,(A) + "M^^M), 



^i(A) 



^?(A) 



(95) 



where 

^<j)p<f>{A) 



2 

TT 



A p1 

dq q(j)Q (q) dx 0o {xq) p {x 
Jo 

1 



W00,2(A) / (ia;cos(solnx)p(x) +W^^^2(A) / (ix sin(so In x) p(x) 
Jo Jo 



(96a) 
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and, similarly, 



^>2(A) 



TT 



dxp{x) + / cos(so Inx)p(x) 
Jo 



<p2(A) f' 



71 



dx sin(so lnx)p{x) 



A A 



n,{A)^^ J J dqdq'p{q,q') 
Here we have used the identity 



A 



2 />! 



c/x p(x) 



0o(a;g) = 4>o{q) cos(solnx) + </)o(g) sin(solnx). 



(96b) 
(96c) 

(97) 



We define those integrals including p{x) as 

dxp{x) = — I dxxe~^^^; 



IT 



9') = — I (ixCOs(So Inx)p(x) — — ^ /^-r nr.cl'e„ in o-W^j PI^- 



TT 



TT 



(ix cos(so lnx)xe 



0^ = — I (ix Sin(So Inx)p(x) — — I cln('e„ in t-W^j 



TT 



TT 



(ix sin(solnx)xe 



(98a) 
(98b) 
(98c) 



The ^j's are constant numbers for a given renormalization condition, and are cutoff inde- 
pendent. Once the parameters a and /3 are numerically fitted, we can obtain the 6'i's from 



Eq. (98). 



Therefore, we now rewrite Sp as 



= 2e,u^^,,{A) + 2e,u^-^^,{K) - 2^^f^ [{e, + e,)<P,{K) + ^3^2(A)] + ^'"^'^^^ 



^i(A) 



^2(A) vr 



(99) 



and see immediately that Sp ~ A^. As illustrated in Subsection VB, p(q',p)'s analytic 



expression (71) can differ from its numerical values ii q/p <t^. 1. One might think that this 



could affect our obtained values of 6'j's since the integrals (98) are from q/p = to q/p = 1 



But p{q,p) is dominated by the region q/p ~ 1, so integrations in the region q/p ^ 1 do 
not affect the results significantly. We have numerically verified that the uncertainty in 
determining 6i is < 1%. 

As the defined summations, Si/j, E2D and Sp, are all ~ A^, we must add them all together 
to cancel the dominant cutoff dependence ~ A^. Therefore, /i2o must satisfy 
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SlD + ^2D + Sp + -/i2o(A)^?(A) = 0, 



(100) 
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FIG. 14: Leading term of the N^LO three-body force, /i2o(A): The sohd (red) hne is the analytic 
expression, and dots are the numerical result. Results are renormalized to reproduce 'jUad = 1-5. 

which yields the result that the leading term of H2, A^/;,2o, is — as expected — of order A^. 
The analytic expression for /120 obtained from Eq. (100) is determined solely by the LO 
renormalization condition, and is thus independent of the additional input required at N^LO. 

We compare this prediction with a numerical calculation of /i2o in Fig- [l4j After repro- 
ducing the atom-dimer scattering length 'yaad = 1-5 at leading order, we have LO parameters 



= —2.012, a = —1.665 and (3 = 1.322, as obtained in Sec. VB for this renormalization 
condition. This yields values of the ^^'s, 61 = -0.038, 62 = -0.036, and ^3 = 0.0088. 
Since the ^j's are small, Sp is generally two orders of magnitude smaller than other terms 



in Eq. (100). Therefore, the 1% error in the determination of ^j's discussed above results in 
only a 10"'^ error in the determination of /i2o- In the numerical calculation of /120; "we also 
fix the NLO and N^LO aad = 1-5/7 ^-nd the N^LO trimer excited-state binding energy shift 
B2 = 0, although the result for /i2o is not affected by the N^LO renormalization condition. 



The numerical values of h2o shown in Fig. 14 agree with the analytic calculation from 



Eq. (100) within an accuracy of 1% for A > 10^7 (the accuracy improves to 0.1% in the 
regime 3 x 10^7 < A < 10^7). This proves that our expression for /120 is an accurate 



prediction. Fig. 14 also shows discrepancies between the analytic and numerical values of 



/i2o at small A. These are due to higher-order 7/A corrections. 

36 



C. The sub-leading cutoff dependence and /i22 



With a fixed 7, /i22 must then cancel the sub-leading cutoff dependence, which is expected 



to be ~ In A in Eq. (39). Since this sub-leading cutoff dependence is two orders lower in 
powers of A compared to the leading one, to extract it from a calculation requires high 
accuracy. In order to obtain a complete analytic expression for /1225 the asymptotic behavior 
of the K- matrix, Koi^q, p; E), must be expanded up to N^LO in all regions {q^ q ^ etc.). 
We are only able to obtain an accurate result for the half-on-shell K-matrix, KQ{k,q; E); 
the near-diagonal {q ~ p) part of the fully-off-shell KQ{q,p]E) was only obtained in an 
approximate form. This approximation likely lacks the high accuracy required for a complete 
calculation of /i22- 

Therefore in this subsection, we first demonstrate numerically that is necessary for 
renormalization. We then perform a partial analytic calculation of /122: we drop the inser- 



tions of NLO terms in Eq. (39), and consider only contributions from the N^LO dimer and 
the N^LO three-body force. This incomplete picture should not infiuence our understanding 
of i?2's sub-leading structure on a power-counting level. In other words, if the sub-leading 
piece of H2 is required to be ~ mE so as to cancel divergences from the insertion of the 
N^LO dimer, we would not expect this feature to be changed by the addition of contributions 
in classes C-E to the calculation. 

Turning first to the numerical calculation, we employ the same parameters as in the 
previous subsection: the atom-dimer scattering length 'jaad = 1-5 is fixed at leading order, 



and we demand that this value is not altered at NLO or N^LO. In Fig. 15 the LO value of the 
trimer excited-state binding energy, B^^\ together with the NLO and N^LO shifts b[^^ and 
are shown as a function of the cutoff A. While B^^^^ and B^'^ are cutoff independent, 
B^^ has a noticeable cutoff dependence, if only one renormalization condition is used, i.e. if 
H2 is assumed to be independent of energy. This demonstrates the necessity for the inclusion 
of an additional counterterm at N^LO. The absence of such cutoff dependence in B^^ shows 
that no energy- dependent counterterm is needed at NLO. 



With the inclusion of the linear- in-energy piece of H2 (73) in the numerical calculation, 
we need an additional parameter at N^LO. Here we fix that parameter by requiring a zero 
N^LO shift of the trimer excited-state binding energy. We then calculate the LO value of 
the trimer ground-state energy, and its NLO and N^LO shifts. These results are plotted 



37 



1.5 



^ 0.5 



PQ 



-0.5 

-1 



-| 1 1 — I — I I I I I 1 1 1 — I — I I I I I 



-| 1 1 — I — I I I I 



10^ 



' □□□□□□pdDDQQDI 



• ^0 

o 5. 



_l I I I ■ ■ ■ ■ I 



_l I I I I I ■ ■ I 



_l I I I I I ■ I 



10 10 
A[Y] 



10 



FIG. 15: The trimer excited-state binding energy as a function of the cutoff A: The dots (dark 
green) are the LO value, the squares (red) are the NLO shifts. The diamonds (blue) are the N^LO 
shifts obtained from using an energy-independent H2. ^Uad = 1-5 is reproduced at all orders. 



as a function of the cutoff A in Fig. [I6| At cutoffs A < 5OO7 cutoff dependence is seen in 
both the NLO and N^LO shift. This cutoff dependence is present because, while such a A 
obeys A ^ 7, A is not well above the typical bound state momentum k. The fact that we 
are examining the deepest bound state means that higher cutoffs must be considered, and, 
indeed, results at all orders converge to a fixed result provided that A ^ 7, /c. We note 
m passmg that the N^LO shift is very large in the "natural units" for this problem, 7^-8^. 
Presumably this is because the presence of the counterterm ~ k means that the N^LO shift 
is really proportional to fc^, not 7^, and is concomitantly larger for this state where k"^ is 
considerably bigger than 7^. We shall return to this issue in the next section. 

To try and better understand the need for a linear-in-energy piece of H2 we define a 
partially summed N^LO K-matrix as 



Kr\k,k;E) 



1 

2^ 



dq . ^^KQ[k,q-E) + —C[k). (101) 



-7 



+ v/3gV4 -mE 



As the 7A dependence is absorbed in the fine tuning of /i20; we consider only the 



0{mE/q^) term in Ko{k,q; E)'s asymptotic expansion (62). Inserting this into Eq. (101), 
we study the cutoff dependence of this new K2{k, k) in an expansion in powers of A, and 
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FIG. 16: The trimer ground-state binding energy as a function of the cutoff A: The dots (dark 
green) are the LO value, the squares (red) are the NLO shifts. The diamonds (blue) are the N^LO 
shifts obtained from using an energy-dependent H2. jaad = 1-5 is reproduced at all orders, the 
trimer excited-state energy's N^LO shift is fixed to zero. 



find: 



+al{k)mE 



1 / X 128/i2o^ 3, \ 32/122^2' 



• (102) 



In order to cancel the leading cutoff dependence (~ A^) of this partial summation, /i2o in 
K^^^^ is required to satisfy: 



x/3 



32/i 



20^2 



$i = 0. 



(103) 



Similarly, /i22 is required to cancel the sub-leading cutoff dependence (~ InA and 0{1)). 



Combining with Eq. (103), we obtain 
1 



2V3 



V3 <Pi \ 2 ' 



(104) 
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FIG. 17: The subleading term of the N'^LO three-body force /122(A): The sohd (red) hne is the 
analytic expression, the dots are the numerical result, ^aad = 1-5 is reproduced at LO, and 



Eq. (106) determines the N LO renormalization condition. 



Therefore, the analytic formula for /i22 in this calculation of the partially-summed K2 is: 

/^22(A) = ^(l + .g ' 



128 ' sin^ [so ln(A /A) - tan'^ sq] 

X {[1 - 3|Ci| cos(argCi)] ln(A/A,) 

-— (sin [2soln(A/A)] -3|Ci|sin [2so ln(A/A) + arg C^] ) 
2so 

1 cos [2so ln(A/A) - tan^^ sq] 



sin [soln(A/A) — tan sq] 



sin [sq ln(A/A) + tan ^ sq] + ^^^^-^ sin [sq ln(A/A) + arg Ci -|- tan ""^So]^! 



X 

\ ^ " ^ ' ' 2 L - ^ . ' " - 

(105) 

where A^, determines the fine tuning of /;,22(A), that is introduced to represent the infrared 
regularization of IA^^^q and W^^^o- •^■i- is independent of the cutoff, but dependent on the 
renormalization condition at N^LO. 



We compare the analytic expression (104) with the numerical result of obtained from 



Eq. (101) in Fig. 17 We again renormalize to 7aad = 1-5 at LO. This determines the value 



A in Eq. (104). Since the partially summed K2 excludes two insertions of NLO terms, it 
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cannot rigorously represent the N^LO corrections to any physical observable. We therefore 
choose two artificial renormalization conditions to determine the unknown constants /i2o and 



part. 



which appear in Eq. (101) for K. 

kr' ^; o) = kr' (0, 0; -B,) = 0. (loe) 

Choosing different N^LO renormalization conditions would yield a different value of the 
parameter in Eq. ( 104). Using numerical results obtained under the N^LO renormalization 
condition (106) we find that the best fit to /i22(A) is obtained with fi^: = 1.752 x 10~^7. The 
remaining behavior of /i22 is a prediction — at least within the context of this subsection's 
approximation — and that prediction agrees with numerical results for /i22 obtained within 
the same truncated calculation of the N^LO atom-dimer amplitude, i^^g"^"^*; accuracy 
of 10-3. 

In this section we have demonstrated the need to include an energy-dependent piece in 
H2 when performing N^LO calculations in the three-body system in SREFT. This was done 
numerically, in a full N^LO calculation, and analytically, in an N^LO calculation that only 
included diagrams of classes A and B. This means that only after reproducing two three- 
body observables in SREFT can we predict other observables at N^LO accuracy. This finding 
agrees with the results of Bedaque et al. in |20j and Barford and Birse in Ref. [25]. However, 
it contradicts the conclusion reached by Platter and Phillips in Ref. [26] . The problem with 
the latter analysis is that the N^LO piece of the dimer propagator was included in the 
modified STM equation, and thereby iterated to all orders. As discussed in Ref. p6], this 
softens the ultra-violet behavior of the atom-dimer K-matrix, i.e. the -ft'(0, g; — -B^) found in 
Ref. falls off more steeply at large q than does ii'o(0,g; —B^)- The integrals computed 
there thus all converge in the limit A — t- 00. A definite A — t- 00 result is found, and no 
additional renormalization appears necessary at N^LO. 

However, as emphasized in Ref. [1^, the existence of a definite A — )■ 00 limit in an EFT 
calculation does not guarantee that the results found in that limit are rigorous consequences 
of the EFT. In particular, in contrast to what was seen in the work of Platter and Phillips, 
the perturbative treatment of ~ rg corrections performed here shows that there is no mod- 
ification of the fall-off at large momenta of the atom-dimer amplitude as that amplitude is 
computed to higher orders in SREFT. The fall off seen at LO (~ 1/g) then determines the 
behavior of integrands in loop integrals to all orders in perturbation theory. We conclude 
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that the approach taken by Platter and Philhps amounts to using the physics at scale tq 
to regulate the loop integrals that constitute the perturbations to the LO result. While the 
results of Ref. [2S] suggest that this is phenomenologically efficacious, it does not constitute 
an EFT calculation if tq ~ /, since in that case other short-distance effects can enter the com- 
putation of atom-dimer scattering at the same scale. Platter and Phillips' non-perturbative 
treatment of the tq corrections could perhaps be justified from an EFT perspective if ro ^ i. 

VII. HELIUM TRIMERS 

To test whether this disagreement between Platter and Phillips' [26] result and our find- 
ing regarding the necessity of H2 at N^LO causes differences in the prediction of three-body 
observables at N^LO, we here apply our full N^LO analysis to observables in trimers of 
Helium-4 atoms. We compare our results with the partially resummed SREFT calcula- 
tions by Platter and Phillips [26] and calculations from realistic potentials by Roudnev and 
Yakovlev [271 EH] and Kolganova et al. [42J. 

It was first observed by Luo et al. [13] that two '^He atoms can form a shallow bound state, 
a dimer. The ^He dimer's scattering length, a = 104l^g A, and binding energy, 3,^ = I.II02 
mK, has been evaluated from measurements of the dimer bond length by Grisenti et al. [S]- 
The effective-range expansion in Eq. ^ can then be used to infer an effective range of the 
order of 10 A. 

However, experimental data for systems of three ''He atoms are still limited. Although 
a three-'^He-atom bound state, i.e. a trimer, has been observed [H], the trimer binding 
energies and the atom-dimer scattering length have not yet been measured. Meanwhile, 
many independent calculations based on realistic potential models [2Z1 [2H1 1121 SSI US] have 
studied ^He trimers. Fairly good agreement has been achieved for trimer binding energies 
and the atom-dimer scattering length. For a review of these calculations, see Ref. [17]. 
These realistic potentials differ at short distances, but all mimic a long-range van der Waals 
potential ~ Cq/t^ for ^He dimers. In Sec.[ll]we used the form of the effective-range expansion 
derived by Gao [31] for such a potential to show that SREFT at N^LO can be used to 
describe systems interacting via such a potential, and only parameters a and tq are needed 
if observables are to be calculated up to an accuracy (9(rQ/a^) or 0{k^rQ). This provides 
an opportunity for a benchmark comparison between calculations of atom-dimer scattering 
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using realistic potentials and SREFT. Such benchmarks have proven very instructive in the 
study of three-nucleon scattering Therefore, in this section we compare our EFT 

results for ^He trimers and the atom-dimer scattering length at 0{rQ/o?) with calculations 
from the TTY potential [271 EH]. 

The two-body parameters a and tq, which are inputs in the SREFT analysis, must then 
be obtained from this potential model. Motovilov et al. found that the TTY potential [51] 
predicts a ^He-^He scattering length of a = 100.01 A and a ^He dimer binding energy 
Bd = 1.30962 mK [13]. From the effective-range expansion, we obtain the atom-atom 
effective range as ro = 7.50(5) A, which provides our EFT expansion parameter ^yr^ = 0.0781. 
We therefore estimate that typical NLO corrections will be ~ 8% and N^LO corrections 
< 1%. 

The SREFT calculations also require one (at LO and NLO) or two (at N^LO) three-body 
inputs. Three-body calculations by Roudnev and Yakovlev from the TTY potential [271 128] 
showed ^He trimer binding energies in the ground and first-excited state of B^^^ = 96.335^ 
and B^^^ = 1.7385^ respectively. The atom-dimer scattering length was obtained as aad = 
1.2057^^. After some initial disagreement, revised numbers for B^^^ and aad from Kolganova 
et al. [^ agree with these values |47j . In what follows we take the results of Ref . [271 EH] , 
to the precision quoted here, as "the TTY results" . 

Upon inserting the N^LO dimer propagator in the modified STM equation. Platter and 
Phillips obtained 5^ ^ = 1.7375(5)5^ and aad = 1.204(1)7"!, where errors in the brackets 
indicate differences from two renormalization schemes. In their calculation, where only one 
three-body input is needed, they can reproduce either B^^^ or aad at LO, NLO and N^LO, 
and predict the other at the corresponding order. Their results are consistent with the TTY 
results with a remaining error ~ 0.1%, consistent with the corrections above N^LO. Platter 
and Phillips also predicted the ^He trimer ground-state binding energy as 5f ^ = 89.45(7)Srf, 
which differs from TTY number by 8%. Since k for this state is large this discrepancy can 
be explained by corrections beyond 0{k'^rl). 

In contrast to these results, we argue that a correct treatment of 0{rQ) effects as perturba- 
tions means that two input pieces of three-body data are needed for proper renormalization. 



In Table. VII , we show our results for ^He trimer observables in two renormalization schemes: 
we fix either or aad to the TTY results at LO and NLO, and use the other as the addi- 
tional parameter at N^LO. If we reproduce aad at all orders, and, in addition, B^^^ at N^LO, 
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we predict the effective range of atom-dimer scattering to be: 



r„^7 = 0.835 + 0.070 + 0.008 = 0.913(1), (107) 

where we have separated the LO, NLO, and N^LO contributions to Vad- This shows the 
excellent convergence pattern for this observable, which is very consistent with the predicted 
expansion parameter rg/a ~ 0.1. Based on this convergence pattern, we conservatively 
predict an N^LO-and-beyond shift in of at most 0.00l7^^. That this is a conservative 
estimate is supported by the difference between the two results for Vad obtained in the 
two different renormalization schemes (fix B^^'^ /aad at LO and NLO), which is an order of 
magnitude smaller: the difference between predictions for rad is 5% at LO, decreases to 0.9% 
at NLO, and to 0.02% at N^LO. 

A different situation is presented by the ^He trimer's ground-state energy. We find (in 



the same renormalization scheme as was used to obtain (107)): 

5f V^d = 97.1 - 7.40 + 27.2 = 116(11). (108) 

Since corrections of the form fc^rg, which are the largest effects for Bf^\ only enter at 
even orders, this may explain the increase in the size of the N^LO correction, compared to 
the NLO one. Here again, the quoted uncertainty on our final result for B^'' comes from 
examining the convergence pattern; while the difference between renormalization schemes 
does not decrease from NLO to N^LO it is markedly smaller than the ~ 8% uncertainty 
we quote here. This uncertainty encompasses the Bf'^ result from the full TTY calculation 
only at the 2a level. Platter and Phillips made a similar argument regarding the accuracy of 
their N^LO result, Bf'^ = 90i?2- In their case, the convergence pattern is less peculiar, and 
the final result obtained reproduces the full TTY calculation within the expected accuracy. 



Eq. (108) suggests that, at best, the SREFT expansion for the trimer ground state is 
poorly convergent. The internal momentum of the state is near the convergence boundary 
of the EFT expansion, making it a questionable venue for the application of SREFT. In 
contrast, a low-momentum observable like Vad seems to converge very well. In order to 
further demonstrate the usefulness of SREFT for low-energy atom-dimer interactions we 
calculated the '^He atom-dimer scattering phase shifts. The results for k cot 6 shown in 
Fig. 18 reproduce aad at all orders, and B^.^^ at N^LO. They should be compared to Fig. 
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where B^^^ is fixed at all orders, and aad at N^LO. In either renormalization scheme (Fig. 
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TABLE L EFT predictions for the "^He trimer binding energies and atom-dimer scattering length 
aad and effective range rad up to N'^LO. Energies and lengths are in units of the dimer binding energy 
and binding momentum respectively. Text in red denotes renormalization conditions (inputs) at 
each order. Results are compared with calculations from the TTY potential \n\ I28|. 
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FIG. 18: ^He atom-dimer scattering phase shifts: aad = 1.2057^^ is fixed at LO (dotted line), and 
also NLO (solid line). Both aad and = 1.738Bd are fixed at N^LO (dashed line). A; cot 5 is in 
units of the ^He dimer binding momentum. 
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FIG. 19: ^He atom-dimer scattering phase shifts: B^^^ = l.lSSBd is fixed at LO (dotted hne), and 
also NLO (sohd line). Both B^^^ and aad = 1.2057~^ are fixed at N^LO (dashed line). kcot6 is in 
units of the ^He dimer binding momentum. 



or Fig. 19), the atom-dimer phase shift converges from LO to N LO. The difference of 



the phase shifts in these two renormahzation schemes is shown in Fig. 20 at each order, 
with the absolute value of this difference, |A {kcot6)\, plotted as a function of the on-shell 
momentum k. |A(A;cot5)| is < 2% at LO, and decreases to less than 0.3% at NLO. The 
N^LO |A(fccot5)| stays below 2 x 10~^ for < k < Based on this renormalization- 
scheme difference, and the convergence pattern of the phase shifts themselves, we deduce 
that N^LO effects from the I = atom-atom amplitude will alter the phase shifts by < 0.1% 
throughout the kinematic range shown. The first atom-atom partial wave with / 7^ that 
is allowed by bosonic symmetry is Z = 2, and this amplitude does not affect the atom-dimer 
scattering amplitude until N^LO. Therefore the phase-shift predictions shown here should 
be accurate to better than 0.2%. 

Refs. |28l SS] give results of atom-dimer phase-shift calculations for the TTY potential. 
However, calculations by Motovilov et al. [45] correspond to an incorrect value of the atom- 
dimer scattering length, aadl = 1-362, a value which was amended to aadl = 1.200(5) 
in later work by the same group of authors [42j. Platter and Phillips used the earlier 
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FIG. 20: Difference between ^He atom-dimer scattering pfiase sfiifts calculated in two renormal- 



ization schemes, which are illustrated in Figs. [18] and 19 respectively. The difference is calculated 
at LO (red dotted hne), NLO (blue dashed hne) to N^LO ( green solid line). |A(A;cot5)| is in units 
of the ^He dimer binding momentum, 7. 



value, and compared their partially-resummed SREFT calculation to the phase shifts from 
Ref. jl5]. Semi-quantitative agreement was found. We employ the correct input aadl = 1.205 
calculated by Roudnev [28j, however, our predicted phase shifts differ from results in Ref. [28] 
by about 3%. Since contributions from higher-partial- wave atom-atom interactions only 
enter at N^LO [M], the origin of this discrepancy is worthy of further investigation. 

High-precision computations with gaussian soft-core atom-atom potentials have recently 
been performed both below |52] and above [53] the three-atom threshold. The soft-core 
potential in these calculations is constructed to reproduce solely the S-wave parameters in 
^He-^He systems (binding energy and scattering length). Such atom-atom potentials are 
supplemented by an energy-independent three-atom force by the same group of authors 
in Ref. |54j in order to obtain a better description of three-atom systems. Since the soft- 
core potential includes only short-distance interactions, calculations with such potentials are 
possiblily equivalent to SREFT calculations. Detailed comparison of SREFT results and the 
calculations of Refs. [52] - [5^ could therefore be quite illuminating. 
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VIII. CONCLUSION 



We calculated the atom-dimer scattering amplitude up to N^LO in the SREFT framework, 
i.e. in a perturbative expansion in r^/a. This amplitude determines three-body observables, 
such as three-body binding energies and atom-dimer phase shifts. We showed the asymptotic 
behavior of LO K-matrices for half-on-shell and fully-off-shell cases. This allowed us to infer 
that, in order to achieve proper renormalization at N^LO, the three-body force at that 
order must contain an energy-dependent part — an inference confirmed by explicit numerical 
calculation. This means that N^LO SREFT calculations in the three-body system require a 
second three-body input, a finding that disagrees with the results of a partially-resummed 
N^LO calculation by Platter and Phillips [26J, where the theory was not renormalized in a 
manner consistent with the perturbative expansion in ro/a and r^k. 

Our perturbative N^LO calculation can be straightforwardly extended to few-body sys- 
tems in nuclear physics, such as the triton or s-wave two-neutron halo nuclei [55| 156] , once 
spin and isospin degrees of freedom are included. For example, the partially-resummed 
N^LO calculation of Ref . [2n] was extended to the three-nucleon system in Ref . [SZ] . Ref . [STj 
found a good result for the triton binding energy, but this was obtained without the energy- 
dependent three-nucleon contact term whose inclusion we are advocating here. Bedaque et 
al. studied three-nucleon observables in SREFT at N^LO, and — motivated by the power 
counting developed in Ref. |TSj — included an energy-dependent three-nucleon contact term 
which played the same role as the piece of the atom-dimer contact interaction [20]. The 
calculation of Ref. [20] also involved partial resummation, which they argued was equivalent 
to a perturbative calculation, up to N^LO accuracy, for momentum-space cutoffs ~ ^/tq. 
However, this partial resummation meant that they could not fully demonstrate proper 
renormalization of their SREFT result at N^LO. Indeed, the only rigorous way to demon- 
strate correct renormalization of the three-nucleon problem is to follow an analysis similar 
to that presented here for three-boson systems. Given the formal similarities between the 
integral equations for three-nucleon system and that for three bosons, we anticipate that the 
conclusions of Ref. pOj (see also Ref. [25]) regarding the necessity of an energy-dependent 
three-nucleon interaction for renormalization at N^LO will be affirmed, once a strict, per- 
turbative N^LO analysis is carried out. 

We have carried out such an analysis for a system of three Helium-4 atoms, calculating 
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trimer binding energies and atom-dimer scattering parameters (flad, f'ad and k cot S) up to 
N^LO and comparing with results from realistic TTY potentials [2H] and the partially- 
resummed EFT result of Ref. [2S]. The convergence of the SREFT expansion is poor for the 
ground-state Helium-4 trimer, since the typical momentum in that state is such that kro ~ 1. 
In the three-nucleon system the triton plays a role analogous to that of this ground-state 
Helium-4 trimer. Thus, slow convergence for the triton in a true N^LO SREFT computation 
(c.f. Ref. [57]) might be expected, because of short-distance effects that scale with fc^rg. In 
contrast, the SREFT expansion converges rapidly for low-energy atom-dimer scattering, 
which facilitates very accurate predictions for the atom-dimer phase shifts in computations 
with the TTY potential. This establishes definitive N^LO predictions for these phase shifts, 
superseding the results of Ref. [26] . Once improved calculations for TTY-potential Helium- 
4-atom-dimer phase shifts become available it will be interesting to compare them to the 
full N^LO SREFT results computed here. 
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Appendix A: Asymptotic expansion of half-on-shell Ko(k,p) 



The half-on-shell K-matrix Ko{k,p; E) obeys Eq. (15) with the on-shell-incoming mo- 
mentum k, off-shell-outgoing momentum p and the 3-body energy E = 3A;^/(4m) — '-f'^/m. 
KQ{k,p; E) is cutoff independent after renormalization. Its result is thus unchanged by use 
of a cutoff A or a different one A' > A: 

Ko{k,p;E) = M{k,p;E) + - [ dq f M{q,p; E)ko{k,q; E) 

7^ Jo -7 + v3gV4 - mE 

= Mik,p;E) + - tdq f M{q,p- E)ko{k,q- E). (Al) 
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By substituting the kernel M{q,p; E) of Eq. (14) into Eq. (Al), we cancel its common part 
on the two sides of the equation: 



A 



2 

IT 

4i^o(A) 
7rA2 



dq 



1 



-7 



In 



mE qp 
q^ko{k,q-E) 



q^ + p'^ + qp — rriE 

g2 _|_ p2 _ _ 







-7 + v/3g74 - mE ttA'^ 



Mk,q; E) 
q^Ko{k,q;E) 







-7+ ^3g2/4-m^ 



(A2) 



As the second integral on the right side of Eq. (A2) vanishes in the limit A' — )■ oo, we have 

,2 1 



In 



-7 + A/3g2/4 - mE qp \q^ + p"^ - qp - mE 



q + p + qp — mE 



Ko{k,q;E) 



4go(A) 
7rA2 



q'Ko{k,q;E) 



(A3) 



-7+ v/3g74-mE' 

Therefore, the effect from the three-body force will cancel the cutoff-dependence in 



Eq. (15), which indicates that 

' k'^ + p'^ + kp 



Koik,p; E) 



kp 



In 



mE 



dq 



/j.2 _|_ p2 _ 



mE 



In 



P -7 + A/3gV4 - mE \q^ + p"^ - qp - mE 



q^ + p^ + qp 



mE 



Koik,p;E) , (A4) 



where the ultraviolet regularization in Eq. (A4) can then be taken to infinity after the 



renormalization. In fact, Eqs. (A3) and (A4) are rigorously valid only if Kq^s asymptotic 



behavior in the region g > A exactly follows Eq. (62). However, this condition cannot be 
obtained numerically using a hard-cutoff regularization. With a finite A, Kq^s asymptotic 
behavior at p ~ A is distorted due to this stiff boundary. We cannot continuously extend 
the asymptotic form of Kq to the region p > A without considering these distortion effects. 

A rigorous way to derive a cutoff-independent integral equation for Kq without introduc- 
ing the three-body force is to apply the subtractive renormalization scheme. For details of 
this method, see Ref. [SH] . 

The function Ko{k,p;E) is bounded when p ~ /c,7; while it is in an expansion of k/p 



and 7/p when p ^ A;,7. This statement is also true for each term in Eq. (A4). 
When p ^ A;, 7, the kernel is expanded as 



k"^ + p"^ + kp — mE 
kp yk'^ + p'^ — kp — mE 



In 



6p^ 



2f 

p4 



+ 



and the propagator is expanded as 
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X 4 mE 
+ — + ■ 



-7 + v^3gV4 - mE ^3 3g 3^3 g^ 3^3 g^ 
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(A5) 



(A6) 



while for arbitrary g, if p 3> k, 7, we have 
+ + qp — mE 



In 



, q^ + p^ + qp 2mE 
In — + 



q/p 



+ 



(A7) 



q^ + p"^ — qp — mE q^ + p"^ — qp p"^ (O'/p)^ + + 1 

In the rest of this section, we simphfy the notation for the half-on-shell K-matrix to 



K{p) = KQ{k,p; E). i.e., its asymptotics at large p is expressed in Eq. (59) as 



k>{p) = p'""-^ + -fDip''°-^ + j^D2p'''>'^ + mE ■ Cip''"- 



(A8) 



which was also derived by Bedaque et al. in pOj. If we compare the inhomogeneous part 



in Eq. (A4) with Ky{p) as regards each term in their 7/p and k/p expansion, we find that 
each term in the inhomogeneous part is 1/p order lower than the corresponding term in 
Ky{p). i.e. corresponds to p^^°~^, and k'^/p'^ corresponds to mE Cip^^°~^, etc.. The 

inhomogeneous terms therefore only affect the overall amplitude of the asymptotic form, 
not the phase. The phase is determined by the factors {D„} and {C„}. Therefore, we drop 
the inhomogeneous part in calculating Kq^s asymptotics, and restore an overall amplitude — 
which is affected by the infrared physics of the inhomogeneous term — in the final calculation. 



After expanding each term in Eq. ( A4 ) , we have 
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111 2 I 2 ^> W 

K<{q) 



dq 
P 



+ ■ ■ ■ 



2mE 



q/p 



-K>{q), 



(A9) 



where /ip separates the integration range into two parts: K^{q) obeys the asymptotic ex- 
pansion at g > fip, and K^{q) is a bounded function at g < fip. 

The bounded function K^{q) guarantees that integrals in the range f^^ are finite, and the 
resultant values of integrals in this region can also be expanded in powers of 7/p and mE/p^. 
At each order in this expansion in powers of the /i-dependence of these parts must 
be canceled by /i-dependence in the high-momentum integral at the corresponding order, 
because is arbitrarily chosen to separate the integral Jq°° in two parts. After canceling the 
/i-dependence, the combination of terms from low- and high-momentum integrations at a 



given order will reproduce the corresponding term in Eq. (A8). 
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Therefore, we can find a series of functions fisg-n{(l), 



fisQ—niSl) 



bounded, if q < ftp 



isQ— n— 1 



if q> HP 



(AlO) 



each one of which matches a term in Eq. (A8). fis^-n a.t q < fip is built to absorb the exphcit 



/i-dependence in Eq. (A9). Therefore, we can replace Ky{q) by /iso-n(Q') in integrals in the 
high-momentum region, and extend the lower limit to 0. After performing the integrals the 



result yields an expression for K^, which can be matched, term-by-term, to Eq. (A8). 



The coefficients, Dn and C„ can be generated from these integrals, which are related to 
Mellin transforms. For example: 

/(iso - n) 



2_Y 



VsJ ULil^ - H^so - k)] 



p{iso-2) + L{iso) 



1 - I{tSo - 2) 

I{s) and L{s) are two types of Mellin transform, that are defined respectively as 

4 



I(s) = 



L(s) 



V3: 



TT Jo 



dx In -—^ — fs(x) 
x^-x + 1'' ^ ' 



V3 



IT Jo 



dx- 



X 



x'^ + x'^ + 1 



9s{x) 



(Alia) 
(Allb) 

(Al2a) 
(Al2b) 



where x denotes the momentum ratio, x = q/p. fs{x) in (A12a) and (A12b) is a dimension- 
less version of the function fs{q)'- 



fsix) 



bounded, if x < fi 



(A13) 



X 



s-l 



if X > 



a definition used in the integrals calculated from the Mellin transform [59] . 

In the rest of Appendix |X| we will explain how we obtain Dn and C„ at each order in the 
7/p and E/p"^ expansion. 



1. So in the leading-order calculation 

At leading order, we have 



i.0-1 4 f dq q^ + p^ + qp 
V Stt Jo P q^ + p- qp 
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By substituting q = xp, we obtain 

4 



dx In 



x^ + X + 1 
x'^ — X + 1 



ISO 



X] 



(A15) 



The integral on the right hand side of Eq. (A 15) is calculated in Ref. [59], which results in 



the analytic expression in Eq. (61). Therefore Eq. (A15) determines the value of Sq by 



I{iso) = 1. 



(A16) 



2. Da's 



By extracting terms proportional to 7 on both sides of Eq. (A4) (and dropping the 
inhomogeneous terms, as previously discussed), we arrive at 



/ 2 \ 4 ^ x^ + x + 1 ^ , , 

Therefore, we relate Di to the Mellin transform I{isQ — 1) by 

2 litso - 1) 



^3 1 - /(zso - 1) ' 



Similarly, we find 



72^2 = 7M ^2 + + 77 ^(^^0 - 2). 



Therefore, D2 is expressed as 
This generalizes to [59] : 
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(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



By inserting the expression for -D„_i back into Eq. (A21) for D„, we find the recursion 
relation 



I{iso - n) 



y/3l{iso - n + 1) [1 - I{iso - n)] 
Using induction, we then derive an expression for Dn. 



Dn-l. 



2 \ " I{isQ — n) 



.VsJ ULA^ - H^so - k)y 

where we have chosen the normalization Dq = 1. 



(A22) 



(A23) 
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3. Ci at 0{mE) 



The coefficient Ci at 0{mE) is derived from 
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By substituting x = q/p, we find that the resuhing Ci is expressed as 



;/(zso-2) + L(zso) 



(A24) 



(A25) 



l-Iitso-2) 

The Melhn transfrom L[s) is computed in Ref. [59], with the final analytic expression given 



in Eq. (61 ). 



Appendix B: Relevant Integrals 



In this section, we define the integrals that are used in the derivation of analytic expres- 
sions for the N^LO three-body forces /i2o and /i22- 

Integrals that contain a single function are defined by 
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(B4) 



where the subscript denotes the power of A present in the integrals for ^„ and n = 1,2. 
Meanwhile, integrals that involve a product of functions are: 
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Lastly, four double integrals used in Sec. jVI 

I - A^ 

^m,i+M) = -J dq<Po{q)^i{q) 



are 
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